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The detection of exotic cosmic rays due to pair annihilation of dark matter particles in the Milky Way halo is a 
viable techniques to search for supersymmetric dark matter candidates. The study of the spectrum of gamma-rays, 
antiprotons and positrons offers good possibilities to perform this search in a significant portion of the Minimal 
Supersymmetric Standard Model parameter space. In particular the EGRET team have seen a convincing signal 
for a strong excess of emission from the Galactic center that has no simple explanation with standard processes. 
We will review the limits achievable with the experiment GLAST taking into accounts the LEP results and we 
will compare this method with the antiproton and positrons experiments. 


The EGRET team fX] have seen a convincing 
signal for a strong excess of emission from the 
galactic center, with 1(E) x E 2 peaking at ^2 
GeV, and in an error circle of 0.2 degree radius 
including the position 1 = 0 ° and b = 0 °. 

This is a particular aspect of a more general 
problem of the diffuse Galactic gamma-ray emis¬ 
sion 2 that also outside the galactic center re¬ 
veals a spectrum which is harder than expected. 
As can be seen in figure ^ the spectrum observed 
with EGRET below 1 GeV is in agreement with 
the assumption that the cosmic ray spectra and 
the electron-to-proton ratio observed locally are 
uniform, however, the spectrum above 1 GeV, 
where the emission is supposedly dominated by 
7r°-decay, is harder than that derived from the 
local cosmic ray proton spectrum [3j. Many dif¬ 
ferent approaches are trying to solve this prob¬ 
lem; one option is the relaxation of the assump¬ 
tion that the local cosmic ray electron spectrum 
is not representative for the entire Galaxy and 
it is on average harder than that measured lo¬ 
cally. Another possibility is that there is some 
variability in the spectral indices of SNR cosmic 
ray sources (for a discussion see TJ). Here we will 
connect the problem of the GeV excess with the 
problem of the missing dark matter in the Uni¬ 
verse and we will examine the possibility to dis¬ 
entangle this effect with the future space 7 -ray 



Figure 1 . Gamma-ray energy spectrum of the 
inner galaxy (300° > l < 30°^) compared with what 
is expected for standard propagation models 
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Figure 2. Gamma-ray energy spectrum of the 
inner galaxy (300° > l < 30°^) compared with what 
is expected for standard propagation models 


and cosmic ray experiments. 

Over the last years our knowledge of the matter 
and energy content of the universe has improved 
dramatically. Astrophysical measurements from 
several experiments are now converging and a 
standard cosmological model is emerging. The 
most significant new data come from recent mea¬ 
surements of the cosmic microwave background 
radiation (CMBR)0 and measurements of the 
Hubble flow using distant supernovae ;6_. 

Current data favor (see for example 7]) a flat 
universe with a cosmological constant = 1 — 
and a total matter density of about 40%=bl0% 
of the critical density of the Universe, with a con¬ 
tribution of the baryonic dark matter less then 
5%. The remaining matter should be composed 
of some yet-undiscovered matter form, such as 
Weakly Interacting Massive Particles (WIMP); a 
good candidate for WIMP’s is the Lightest Super- 
symmetric Particle (LSP) in R-parity conserving 
supersymmetric models. 

The motivation for supersymmetry at an ac¬ 
cessible energy is provided by the gauge hierar¬ 
chy problem [HJ, namely that of understanding 
why mw rap, the only candidate for a fun¬ 


damental mass scale in physics. This difference 
introduces problems because one must fine-tune 
the bare mass parameter so that is almost ex¬ 
actly cancelled by the quantum correction in or¬ 
der to obtain a small physical value of mw • This 
seems unnatural, and the alternative is to intro¬ 
duce new physics at the TeV scale and to postu¬ 
late approximate supersymmetry ||, where boson 
and fermion partners naturally induce cancella¬ 
tions in quantum corrections in case 

| m 2 B — m 2 F | < ITeU 

This is also the reason to expect that, if super- 
symmetry is real, it might be accessible to the cur¬ 
rent generation of accelerators and in the range 
expected for a cold dark matter particle. 

The minimal supersymmetric extension of the 
Standard Model (MSSM) pH) has the same gauge 
interactions as the Standard Model; most often 
its phenomenology is simplified by reducing its 
parameter space to only six parameters: the higgs 
mixing parameters /a that appears in the neu- 
tralino and chargino mass matrices, the common 
mass for scalar fermions at the GUT scale mo, the 
gaugino mass parameter Mi/ 2 , the trilinear scalar 
coupling parameter A, the ratio between the two 
vacuum expectation values of the Higgs fields de¬ 
fined as tan/5 = V 2 /V 1 =< H 2 > / < Hi > and 
the mass of the pseudoscalar Higgs 

The LSP is expected to be stable in R-parity 
conserving versions of the MSSM, and hence 
should be present in the Universe today as a cos¬ 
mological relic from the Big Bang m- R-parity 
is a discrete symmetry related to baryon number, 
lepton number and spin: 

R = (-1) 3£?+L+25 

It is easy to check that R=+l for all Standard 
Model particles and R=-l for all their supersym¬ 
metric partners. There are three important con¬ 
sequences of R conservation: (i) sparticle are al¬ 
ways produces in pairs; (ii) heavier sparticles de¬ 
cay into lighter sparticles and (iii) the LSP is sta¬ 
ble because it has no allowed decay mode. 

The LSP is expected also to be neutral, be¬ 
cause with an electric charge or strong interac¬ 
tion, it would have condensed along with ordinary 
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baryonic matter during the formation of astro- 
physical structures, and should be present in the 
Universe today in anomalous heavy isotopes na- 
This leaves as candidates a sneutrino with spin 0, 
the gravitino with spin 3/2 and the neutralino x 
that is a combination of the partners of the 7 , Z 
and the neutral Higgs particles (spin 1/2). 

Searches for the interactions of relic particles 
with nuclei rule out sneutrinos with masses lighter 
than few TeV l I 2], while the gravitino could con¬ 
stitute warm dark matter with a mass around 1 
keV. So the best candidate for cold dark matter 
appears to be the neutralino y. The experimental 
LEP lower limit on m x is M 

rn x > 50 GeV 

As m x increases, the LSP pair annihilation rate 
in dark structures decreases, but, as we will show 
below, up to ^ 300 GeV, the distortion of the 
spectrum of the diffuse gamma ray background 
due to a neutralino induced component up to the 
neutralino mass can be a possible signature of the 
existence of the LSPHSJ. 

How can this kind of signal be seen? 

In figure [2| is shown the EGRET data from 
the Galactic center, the diffuse gamma ray back¬ 
ground flux expected from the standard interac¬ 
tions and propagation models of cosmic ray pro¬ 
tons and electrons and an example of the flux 
due to neutralino annihilation in the dark mat¬ 
ter halo. In this case the signal is for a ^ 100 
GeV neutralino and for the bb annihilation chan¬ 
nel (the spectral shape of the other channels is 
very similar). 

The total flux is: 

(f)ry = NjjfPl) + N X (f) X , ( 1 ) 

where and N x are normalizations parameters 
for the standard and exotic flux, respectively. N\ 
is associated to the interstellar medium column 
density n(Z), integrated along the line of sight. 

The background flux can be written in the fol¬ 
lowing way: 

< j>b(E 7 ) = Em(E 7 )/ ( cm 2 sr ) ( 2 ) 

N b= f j-n{l)/ (cm _ 2 sr _1 ) (3) 

Jl.o.s 47r 


where Em(E^) [GeV~ 1 s~ 1 ] is the emissivity per 
hydrogen atom, which gives the number of sec¬ 
ondary photons with energy E 1 , emitted per unit 
time and unit energy. The two main production 
mechanisms for the background are through the 
7T° production and the inverse compton scatter¬ 
ing. We have considered two channels for the i r° 
production, i.e. primary p and He. The processes 
are of the type: 

P + X —► .. —> 7T° —> 2y, (4) 

HT-bV—>7T°—>2 7 

where X can be interstellar H or He nucleus. 
We have simulated the spectrum using the Gal- 
prop computer code |T 6 | . The WIMP flux cj) x in 
eq. © originates from the annihililation of a cou¬ 
ple of generic WIMPs, of mass m x , in allowed 
tree-level final states, of which the leading chan¬ 
nels are usually fefe, cc, ££, IV + IV - , Z°Z°. Light 
quark states are suppressed by m 2 /ra 2 in the an¬ 
nihilation cross section, where m/ is the inter¬ 
mediate fermion mass. Higgs bosons annihilation 
states are included, allowing for their decay into 
particles for which we do simulate. The contin¬ 
uum gamma ray flux coming from a direction that 
forms an angle ^ with the direction of the Galac¬ 
tic center can be written as: 

* E ,*) = ^L/wm ( 5 ) 

where m x is the WIMP mass, av is the total anni¬ 
hilation cross section times the relative velocity of 
the two annihilating WIMPs (in the galactic halo 
v/c ~ 10 -3 ) and p(l) is the WIMP mass density 
along the line of sight. Contributions along the 
line of sight are then summed. The equation © 
can be factorized into two pieces, one depending 
only on the cross section and the WIMP mass, 
and the other one depending on the WIMP distri¬ 
bution in the galactic halo. This factor depends in 
a crucial way from the halo model used. We now 
write the differential flux <j) x of eq. ©, in units of 
era - 2 GeV -1 , as: 

= 1.87 •10- 11 f in _ 29 2<n, 3 X 
\10 z ^cm 6 s 1 J 

/ 10 GeV V f dN 
V m x J \ dE 
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where dN/dE is the the number of photons pro¬ 
duced per unit energy in each intermediate an¬ 
nihilation channel. We then introduced the di¬ 
mensionless function J(V ; ) that encodes the de¬ 
pendence from the halo density profile: 

■w-sMo5d^)P<w>< 7 > 

In such a way the adimensional normalization fac¬ 
tor N x of eq. © is then precisely J( ip). Actually 
we have averaged the value of J($) over a solid 
angle AU around the direction determined by the 
angle ip: 

W)> AO = Xo / (8) 

where we have chosen the same EGRET region 
around the galactic center of AQ = 2.15-10' 3 sr. 

As it can be seen from figure O the fit to the 
data greatly improve when a neutralino compo¬ 
nent is added. Of course this cannot be assumed 
as evidence for the existence of supersymmetry 
particles as the dark matter component of the 
halo, but as an indication that more experiments 
with greater sensitivity and exposure are needed. 

In the next session we present one future pos¬ 
sibility, i.e. the experiment GLAST. 

1. The Gamma-ray Large Area Telescope 
GLAST 

The standard techniques for the detection of 
gamma-rays in the pair production regime energy 
range are very different from the X-ray detection. 

For X-rays detection focusing is possible and 
this permits large effective area, excellent energy 
resolution, very low background. For gamma-rays 
no focusing is possible and this means limited 
effective area, moderate energy resolution, high 
background but a wide field of view. This pos¬ 
sibility to have a wide field of view is enhanced 
now, in respect to EGRET, with the use of silicon 
detectors, that allow a further increase of the ra¬ 
tio between height and width, essentially for two 
reasons: a) an increase of the position resolution 
that allow a decrease of the distance between the 
planes of the tracker without affecting the angu¬ 
lar resolution, b) the possibility to use the silicon 


Figure 3. The GLAST instrument, exploded to 
show the detector layers in a tower, the stacking 
of the Csl logs in the calorimeter, and the inte¬ 
gration of the subsystems. 


detectors themselves for the trigger of an events, 
with the elimination of the Time of Flight sys¬ 
tem, that require some height. The Gamma-ray 
Large Area Space Telescope (GLAST) ESI, has 
been selected by NASA as a mission involving 
an international collaboration of particle physics 
and astrophysics communities from the United 
States, Italy, Japan, France and Germany for a 
launch in the first half of 2006. The main sci¬ 
entific objects are the study of all gamma ray 
sources such as blazars, gamma-ray bursts, su¬ 
pernova remnants, pulsars, diffuse radiation, and 
unidentified high-energy sources. Many years of 
refinement has led to the configuration of the ap¬ 
paratus shown in figure 03 where one can see the 
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Figure 4. Total photon spectrum from the galac¬ 
tic center from standard propagation models and 
from one neutralino annihilation models and the 
kind of statistical errors that it is expected in three 
years with GLAST. 


4x4 array of identical towers each formed by: •Si- 
strip Tracker Detectors and converters arranged 
in 18 XY tracking planes for the measurement 
of the photon direction. • Segmented array of 
CsI(Tl) crystals for the measurement the photon 
energy. • Segmented Anticoincidence Detector 
(ACD). The main characteristics, shown in fig¬ 
ures 0 to 0 are an energy range between 20 MeV 
and 300 GeV, a field of view of ^ 3 sr, an energy 
resolution of ^ 5% at 1 GeV, a point source sen¬ 
sitivity of 2xl0 -9 (ph cm -2 s -1 ) at 0.1 GeV, an 
event deadtime of 20 /as and a peak effective area 
of 10000 cm 2 , for a required power of 600 W and 
a payload weight of 3000 Kg. 

The list of the people and the Institution in¬ 
volved in the collaboration together with the 
on-line status of the project is available at 
http://www-glast.stanford.edu A description of 
the apparatus can be found in m and a descrip¬ 
tion of the main physic items can be found in 

EQ- 

In figure 0 is shown the total photon spectrum 
from the galactic center from standard propaga- 



Figure 5. GLAST effective area as a function of 
energy including all background and track quality 
cuts compared with EGRET’s one 



Figure 6. GLAST angular resolution. 



Figure 7. GLAST relative effective area for 1 
Ge V photons as a function of the zenith angle 







6 



Figure 8. GLAST and EGRET energy resolution 


tion models and from one neutralino annihilation 
models and the kind of statistical errors that it 
is expected in three years with GLAST in the 
case of a moderate value of 500) which 

is within the allowed ranges of both the modified 
isothermal and cuspy halos. This effort will be 
complementary to a similar search for neutralinos 
looking with cosmic-ray experiments like the next 
space experiments PAMELA [Ti and AMS^, at 
the distortion of the secondary positron fraction 
and secondary antiproton flux induced by a signal 
from a heavy neutralino. 

2. Conclusion 

The gamma-ray space experiment GLAST is 
under construction. Its time of operation and 
energy range is shown together with the other 
space X-ray satellite and gamma-ray experiments 
in figure fTTl Note that it will cover an interval not 
covered by any other experiments. Note also the 
number of other experiments in other frequencies 
that will allow extensive multifrequency studies. 
In the last decade, ground-based instruments 
have made great progress, both in technical and 
scientific terms. High-energy gamma rays can be 
observed from the ground by experiments that de¬ 
tect the air showers produced in the upper atmo¬ 
sphere. In fTSlthe GLAST sensitivity is compared 



Figure 9. Distortion of the secondary antiproton 
flux induced by a signal from a heavy neutralino. 
The PAMELA expectation in the case of exotic 
contributions are shown by black squares 


PAMELA expectation 
for 3 years of operation 
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Figure 10. Distortion of the secondary positron 
fraction The PAMELA expectation in the case of 
exotic contributions are shown by black squares 
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Energy versus time for X and Gamma ray detectors 



Figure 11. Timeline schedule versus the energy 
range covered by present and future detectors in 
X and gamma-ray astrophysics. 



Figure 12. Sensitivity of present and future de¬ 
tectors in the gamma-ray astrophysics. 


with the others present and future detectors 
in the gamma-ray astrophysics range is shown. 

The predicted sensitivity of a number of opera¬ 
tional and proposed Ground based Cherenkov 
telescopes, CELESTE, STACEE, VERITAS, 
Whipple is for a 50 hour exposure on a single 
source. EGRET, GLAST, MILAGRO, ARGO 
and AGILE sensitivity is shown for one year of all 
sky survey. The diffuse background assumed is 
2 • 10 -5 photons cm~ 2 s~ x sr ~ l (100 MeV/E ) 1A , 
typical of the background seen by EGRET at 
high galactic latitudes. The source differential 
photon number spectrum is assumed to have 
a power law index of -2, typical of many of 
the sources observed by EGRET. Above lGeV 
the GLAST sensitivity shown is not limited 
by background, but rather by the requirement 
that the number of source photons detected is 
at least 5 sigma above the background. The 
AMS sensitivities is from 123 • Note that on 
ground only MILAGRO and ARGO will ob¬ 
serve more than one source simultaneously. The 
Home Pages of the various instruments are at 
http:// www-hfm.mpi-hd.mpg.de/ CosmicRay/Cosmf 
As is shown in inn, indirect dark matter searches 


and traditional particle searches are highly com¬ 
plementary and, in the next few years, many 
experiments will be sensitive to the various po¬ 
tential neutralino annihilation products. These 
include under-ice and underwater neutrino tele¬ 
scopes, atmospheric Cerenkov telescopes, high 
altitude extensive air showers detectors and 
the already described space missions GLAST, 
PAMELA together with AMS. In many cases, 
these experiments will improve current sensitivi¬ 
ties by several orders of magnitude and probably, 
as it is shown in EH,” all models with charginos 
or sleptons lighter than 300 GeV will produce 
observable signals in at least one experiment in 
the cosmologically preferred regions of parameter 
space with 0.1 < Q x h 2 < 0.3 ” before LHC. 
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